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ABSTRACT 

Context. 

Aims. We investigate the accuracy of astrometric measurements with the VLT/FORS1 camera and consider potential applications. 
Methods. The study is based on two-epoch (2000 and 2002/2003) frame series of observations of a selected Galactic Bulge sky 
region that were obtained with FORS 1 during four consecutive nights each. Reductions were carried out with a novel technique that 
eliminates atmospheric image motion and does not require a distinction between targets and reference objects. 
Results. The positional astrometric precision was found to be limited only by the accuracy of the determination of the star photocentre, 
which is typically 200-300 //as per single measurement for bright unsaturated stars 6=18- 19. Several statistical tests have shown 
that at time-scales of \-A nights the residual noise in measured positions is essentially a white noise with no systematic instrumental 
signature and no significant deviation from a Gaussian distribution. Some evidence of a good astrometric quality of the VLT for 
frames separated by two years has also been found. 

Conclusions. Our data show that the VLT with FORS 1/2 cameras can be effectively used for astrometric observations of planetary 
microlensing events and other applications where a high accuracy is required, that is expected to reach 30-40 yuas for a series of 50 
frames (one hours with6 specia i filter). 

Key words, astrometry - atmospheric effects - instrumentation: high angular resolution - planetary systems 



1. Introduction 

The provision of high-precision astrometry (of a few 100 yuas 
and below) is much desired, but so far it is not widely available. 
Over the recent years, photometric microlensing has proven its 
feasibility for not only detecting massive gas giant planets (Bond 
et al. 2004; Udalski et al. 120051 1, but even cool rocky/icy plan- 
ets (Beaulieu et al. 2006). A fourth planet, of roughly Neptune- 
mass, was most recently claimed by Gould et al. (2006). Several 
further potential detections of planets by microlensing have how- 
ever been missed (e.g. Jaroszynski & Paczynski 2002) since a 
proper characterization requires good data quality and a dense 
sampling. In fact, it is necessary to overcome ambiguities and 
degeneracies of binary lenses themselves (Dominik 1999), be- 
tween mass ratio and source size (Gaudi 1997), and between 
binary-lens and binary-source systems (Gaudi 1998), as shown 
explicitly by Gaudi & Han (2004). The observation of the shift 
of the centroid of light composed of the lens star and the images 
of the source star with time, so-called astrometric microlensing 
(H0g et al. [19951 Hard y & W alker [19951 Miyamoto & Yoshii 
1980; Dominik & Sahu 2000), provides a two-dimensional vec- 
tor, whose absolute value is proportional to the angular Einstein 
radius (comprising information about the lens mass and relative 
lens-source parallax), in addition to the magnification provided 
by the photometric observations, where the latter only allows 
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to measure the angular Einstein radius if the finite size of the 
source star can be assessed in the light curve. While Safizadeh et 
al. ([1999) first discussed the possibility of using astrometric mi- 
crolensing for discovering and characterizing extra-solar planets, 
Han & Chunguk ( 2002) found that the planetary deviation to the 
astrometric signal expected from a Galactic bulge star lensed by 
a foreground star and a surrounding Jupiter-mass planet is about 
10-200 /ias. Rather than requiring a stability of the detector on 
time-scales that correspond to the orbital period of planets as for 
the detecting of the astrometric shift seen in its host star due to 
motion around the common center-of-mass (e.g. Quirrenbach et 
al. 12004; Pravdo et al. 2005), the duration of the microlensing 
perturbations is just between a few hours and a few days. 

Since ground-based monopupil telescopes were thought to 
be limited by atmospheric noise and therefore have rather poor 
precision (Lindegren 1980), only the use of high-precision opti- 
cal interferometers such as SIM (Boden et al. [T998 ; Paczynski 
1998), Palomar Testbed Interferometer (Lane & Muterspagh 
|20Q4l l Keck (Colavita et al. [19981 , and VLTI (Wittkowski 
et al. 120041 1 was considered so far for carrying out observa- 
tions with such precision. A quick progress in adaptive optics, 
and its application to astrometry (Pravdo et al. 2005) leads to 
impressive results. However, the theoretical considerations by 
Lazorenko & Lazorenko (120041 1 and a first set of correspond- 
ing data (Lazorenko 2006 hereafter Paper I) demonstrated that 
astrometric measurements with very large ground-based tele- 
scopes are not atmospheric-limited. Instead, a high precision be- 
yond this limit is achievable by elimination of atmospheric im- 
age motion at the data-reduction phase. For that purpose, a field 
of reference stars is configured into a virtual high-pass filter that 
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absorbs most significant low-frequency modes of the image mo- 
tion spectrum. For telescopes with a fully illuminated entrance 
pupil, the atmospheric error <x at rapidly decreases with the in- 
crease of the telescope entrance pupil D. At D — > oo, the asymp- 
totic dependence is 

<T*=B k R n l 6 D-' il2 T- X12 (1) 

where T is the integration time, R is the angular reference field 
radius, and is the filter parameter. This relation involves 
an additional factor ^R/D as compared to that predicted by 
Lindegren ( 119801 ), so that large telescopes with D = 8 - 10 m 
and above are favourable. A specific amplitude apodization to 
the entrance pupil of future telescopes may further suppress the 
atmospheric error to <x at ~ /? fc / 2 £)-*/ 2+1 /3j'-i/2 where k is the re- 
duction parameter. 

The astrometric performance of the VLT on very short time- 
scales could have been estimated from test reductions of a four- 
hour series of FORS2 frames obtained in the region of the 
Galactic bulge (Paper I). A positional accuracy of 300 fias per 
single R-16 mag star measurement at 17 sec exposure was 
shown to be limited only by the number of photons received 
from the observed star. 

The current study extends the certification of high-accuracy 
astrometric performance with the VLT to time scales of 1-4 
days, which is sufficient for monitoring astrometric planetary 
deviations in microlensing events. In Sect. 2 of this paper, we 
lay out the observational basis of this investigation, describe 
the image processing, and show how the stellar centroids are 
computed. Particular problems in astrometric reductions that are 
based on background reference stars with unstable positions, 
and a new reduction method are discussed in the subsequent 
Sect. 3. Reduction of the observational data with this method 
is described in Sect. 4. Sect. 5 focusses on the noise and cor- 
relations in the measurements that, if exist, can be related to 
the different parts of the telescope (optics, mechanics, electron- 
ics), to the atmosphere, etc., preventing statistical improvement 
of the accuracy by means of simple accumulation of the number 
of frames. 

2. Observations and computation of centroids 

For our study, we were able to use VLT frames that are available 
in the ESO Archive: those obtained with the FORS 1 camera of 
UT1 in the sky region near the neutron star RX J0720.4-3 125 un- 
der the ESO program 66.D-0286 (Motch et al. 120031 were found 
to be best-suited for our purposes. Images were acquired during 
4 nights from 20 to 23 December 2000 (40 frames) and 4 nights 
from 29 December 2002 to 02 January 2003 (25 frames). During 
each 5 to 12 frames were obtained in B filter with ~ 10 mn ex- 
posure time and FWHM varying from 0.45 to 0.85", where the 
high-resolution mode with 0.10"/px scale and about 3.3 X 3.3' 
field-of-view was used. 

These data enable us to carry out a detailed and comprehen- 
sive study of the astrometric precision that is achievable with the 
VLT/FORS1 at short time scales. Moreover, our analysis serves 
as a test reduction for future observations. Unfortunately, be- 
tween some frames, the images turned out to be displaced by 
more than 200 px. For that reason, peripheral stars frequently 
needed to be omitted for much jittered frames. Moreover, for as- 
trometric purposes, the choice of the B filter is not optimal, since 
it results in larger residual atmospheric chromatism as compared 
to /^special- However, it were just these unfortunate circumstances 
that forced us to look for a reduction that significantly improves 
the quality. 



Observations were obtained with the Longitudinal 
Atmospheric Dispersion Corrector (LADC), which is a 
two-prism optical device (Avila et al. 1 19971 1 that allows to 
correct for atmospheric differential chromatic refraction (DCR). 
DCR is caused by the atmosphere, which acts like a prism with 
a certain dispersion and thereby smears and displaces stellar 
images (e.g. Monet et al. |1992t . The size of this effect depends 
on the zenith distance and the colour of the observed star and 
can be modelled as shown in Eq. ( TTZb . 

Raw images were debiased and flat-fielded using master cal- 
ibration frames taken from the ESO Archives. Positions x, y of 
stellar centroids were computed using the profile fitting tech- 
nique that had already been applied earlier for the high-precision 
astrometric reduction of FORS2 images (Paper I). Stellar pro- 
files in 10 x 10 px windows were fit by a sum of three modified 
elliptic Gaussians involving 12 free model parameters altogether. 
The dominant Gaussian centered at x, y with extent parameters 
cr x c , cr c along corresponding coordinate axis and futher speci- 
fied by its orientation angle, was considered to contain a flux / 
mounting to about 2/3 of the total light received from the star. 
Other auxiliary Gaussians allowed profile fitting to the photon 
noise limit at which residuals of pixel counts from the model are 
characterized by x 2 ~ L0/px (except for central parts of im- 
ages of bright stars). Objects with^ 2 > 3.0/px deviations were 
rejected and only unsaturated images with B > 18 mag were 
processed. 

Due to the complex shape of the PSF, the precision e of 
the image centroiding was obtained by numerical simulation, 
for which a set of randomly generated star images with profiles 
comparable to the observed ones was created and where Poisson 
noise was added to pixel counts. This simulation shows that s 
can be well-approximated by 

e= — %(1 +a r 0J ), (2) 
2.46 V/ 

where 9 denotes the seeing (which turned out to be 8 = 3.10<x G 
in our three-component model), / refers to the electron count, 
and ao is a characteristic constant that depends on the back- 
ground noise. For the model background fluxes bo of 810 or 
540 e~/px (average values at 2000 and 2002-2003 epochs), ao 
was found to be 1150 or 887, respectively. Eq. becomes an 
exact relation only for a fixed average seeing 6q = 0.608" at 
which ctq — 1 .96 px, but it provides a reliable approximation for 
variable FWHM of that order. 

3. The reduction process 

3.1. Filtration of atmospheric image motion 

Let us shortly recall the previously discussed process of atten- 
uation of the image motion spectrum (Lazorenko & Lazorenko 
2004), which is based on the virtual symmetrization of the ref- 
erence field. Consider some region of the sky with N stars being 
imaged in m = 0, 1 . . . M consecutive frames, where x/ m , y\ m are 
the centroids of / = 1,2... N m stars measured in frame m. For 
various reasons (jittered images, cosmics, outliered data, etc.) a 
significant fraction of centroids cannot be obtained from every 
frame, therefore N m < N. Here and below, indices i, j, I refer to 
the target, reference, and star of any type, respectively, and m - 
is assigned to the reference frame. Upper indexes x, y, if used, 
refer to the coordinate axes. A list of variables most frequently 
used in the Paper with a short description is given in Table Q] 

The differential position Xj m of the 2-th target star in the m-th 
frame is obtained from CCD coordinate differences x> m - x, m of 
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Table 1. A list of frequently used variables 





definition 


first use, 




equation 


x,y 


measured star positions 


(3) 


X, Y 


positions x, y filtered of image motion and given 
in the system of local reference stars with 

fivi^ri nncitinnc 
11ACU U(J A 1 LI Ul IS 


(3) 


A\ X> 


same as X, Y but given in a single all-frame 






system with "floating" reference stars 




e 


displacement (average for a given frame series) 
of the star from its position in reference frame, 
for astrometric reasons (proper motion etc.) 


(14) 


i\i y 


same as £ x , ^ y but given with respect to the 
mathematical expectation of the star position in 
reference irame 


(10) 


p 


chromatic coefficient; atmospheric displacement 
oi tne star image is proportional to p 


(12) 


d 


r* n ivvm n 1 1 r* f*{~\f*rnr*if*nt' npcr*i"i npc r*nmnpncnfinrr 
CI 11 (JlllclLLU L.UC111U1C11L, U-CSUllUCa UUlllUCllaaLlllg 

action of the LADC 


(191 


h 


differential chromatic displacement (DCR effect) 
of the star image 


(12) 


M 


atmospheric image motion effect in x, y 


(10) 


AM 


residual atmospheric image motion in X, Y 


(14) 


e 


error of photocenter measurement 


(10) 


e 


variance of e 


(1,25) 


u 


residuals of conditional equations (23) 


(26) 


N 


noise in observations that enters Eqs.(23) 


(17) 


a- 2 


observed variance of N found from residuals u 


(26) 



the target and field star centroids. More specifically, it is given 
as the weighted average 

— ^ ^ Uij(Xj m Xjm) ~ %im ^ ^ @ijXjm> 111 — 1,2... M (3) 



7=1 



7=1 



of these differences, where the weights ay satisfy the normaliza- 
tion condition 



(4) 



7=1 

and are chosen, so that the effect of image motion effect in X m is 
minimized. The prime indicates that the summation is carried out 
over the subset £!,„, of stars used as reference for z'-th star only, 
where the index j = z is omitted. The corresponding positional 
y-coordinate F„„ is obtained in analogy to Eq. ([3]). 

Atmospheric image motion displaces each position x, m by 
a small angle M m , which results in X m being displaced by 
AMi„, = Yjj' a ij(Mi,„ - Mj m ). The spectral power density G(q) 
of AMi m in the space of spatial frequencies q is the product of 
two factors. The first factor F'(q) depends on D, exposure time 
T, and properties of atmospheric turbulent layers generating im- 
age motion. The second factor depends only on the geometry 
of background star positions relative to the target i, and can be 
expanded in a series of even powers of q. Thus, 



G(q) = F'(q) ^ q 2s F 2s (xi, yr, X\ , yi . . . x N , y N ) , 



(5) 



s=1 



where x, y are positions not affected by image motion (at zero 
turbulence). The coefficients F 2 s are quadratic functions of co- 
ordinate differences that fastly decrease with s and take the form 



(6) 



(2. s) 

where A „ are the characteristic constants. With the choice of 
coefficients fly, all sums in brackets can be turned to zero for 
a specific A ~jj . In this way, the leading term F 2 and, possibly, 
several subsequent F 2s terms are eliminated, so that G(k) will 
depend largely on the first non-zero component of some op- 
tional mode k = 4, 6 For this purpose, the coefficients fly 

should satisfy each of k(k + 2)/8 - 1 conditions 



S'flytao - XjoTiyiQ-yjof 
j 



0, a+j3= 1 



1. 



(7) 



where k > 4, a and ft are non-negative integers, and positions 
x, y are substituted by respective measured positions of stars in 
the reference frame. The coefficients fly are determined for each 
target i by solving Eqs. with normalizing condition (0, pro- 
vided that at least N' = k(k + 2)/8 reference stars are available. 
Because usually N m » N', the system of equations (0 is re- 
dundant and solvable with some useful restrictions on ay, for 
example Eq. ( TT3T >. Due to the elimination of principal modes of 
the image motion spectrum, the resulting values of X lm are char- 
acterized by the residual image motion AAt; m that involves only 
uncompensated spectral modes of orders s > k/2 and therefore 
is much smaller than its initial value. Approximately, the vari- 
ance of residual image motion is expressed by Eq. CD, where R 
is to be substituted by the radius Rj of local reference field for 
z'-th target. The advantage of using larger orders in k is reflected 
in the coefficient which decreases with k. In particular, for 
k — 12, its value is only 20% of that measured at the VLT with 
k = 4 (Paper I). 

Next, let us introduce polynomial base functions fu = 1, 

hi = xio, fu = yio, ■ ■ ■ fwl = *Jyf of coordinates x m , y m 
referring to the /-th star in the reference frame, with indices 
w = 1 , 2 . . . N' equal to the sequential number of combinations 
of a and B in equations (O. Taking advantage of the equality 
Z/dijiXio-XjoTfyio-ypf = xf yf -2/ayx^ , we canrewrite 
equations (|4]i and (0 in terms of functions f w i as 



2'fly/i / = 1, 2'ay/w./ = fwi, w = 2, 3 . . .N' 
j j 



3.2. Conversion from X, Y to displacements A x , A y 



(8) 



The image displacement X m of the z'-th star between frame m 
and the reference frame, filtered from low-order image motion 
and systematic distortions of the field, is measured relative to 
a set Q,'„, of background stars whose positions are assumed to 
be fixed. However, since such are not available, we bypass this 
limitation by introduction of astrometric displacements A that 
are related to all stars in the frame without distinction between 
reference and target objects. We consider X m and A related by 



Xjn, 



XlO + Plm + Al„ 



(9) 



associating the measured position of stars in frame m with their 
position in the reference frame. Here, P/ m = £ w c wm f w i is a poly- 
nomial function that contains powers of x, y not higher than 
k/2 - 1 and measured at discrete set of points xio, yio. Via co- 
efficients c wm , this function describes the systematic variation of 
stellar positions with the change of optical aberration and other 
similar effects. The second quantity 



Ai m =ii + e/m - e /0 + hi,„ - hio + Mi,,, - MlO 
is orthogonal to P; m and satisfies conditions 
Zi *Uwi = 0, Zi = 0, w = 1,2 . . .N' 



(10) 



(11) 
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In fact, Eqs. (0 and ( fTTT i give an expansion of x\ m into power se- 
ries of f w i, with Ai m representing the remainders. A\ m includes the 
"science" astrometric signal (sum of proper motion, parallax, 
probable microlensing and reflex motion caused by the planetary 
companion) measured relative to the reference frame, random er- 
rors of photocenter measurement e/ m , the atmospheric chromatic 
displacement /i/,„ corrected by LADC, and the atmospheric im- 
age motion Atj». These quantities (or their combinations) are 
also orthogonal to f w i, with their low-order w < N' expansion 
terms included in P\ m via coefficients c wm . 

The components of h/ m along the x- (opposite to right ascen- 
sion) and y- (along declination) axes are given by (Paper I) 



for x - axis 
for y - axis 



(12) 



h* = pi tan Zm sin y m + di tan £ m sin y m 

Kn = ~pi tan z »< cos y™ ~ d > tan (»t cos y* 

Here, p\ is the coefficient of differential chromatic displacement 
for the Z-th star, di is a negative displacement produced by LADC 
to compensate for DCR effect, approximately of the same abso- 
lute size but opposite sign, z m is a zenith distance at the midpoint 
of exposure of the m-th frame, y m is the parallactic angle, and £ m 
is the zenith angle of the LADC setting, which remains the same 
for a single series of observations. 

Substitution of Eqs. (O for Xj m and xj m into Eq. (O leads 
to elimination of the polynomials P\ m , while fulfilling the condi- 
tions given by Eq. ([8]). Taking into account that X® = £ j<hfa«y- 
Xjo) = due to Eqs. ©, we obtain 



At 



- ^ 0,ijXj m - Xu 



(13) 



An analogous expression for transformation of Y- m to results 
for measurements along the y-axis. Eqs. ( fT3t . for each frame m, 
form a system of i — 1 . . . N m linear equations with N m unknowns 
Af , The solution for a given frame m is found irrespective of so- 
lutions for other frames. The matrix of the system ( fT3b . however, 
turns out to be degenerated due to N' linear dependences (Q be- 
tween its columns (coefficients fly), and consequently it can be 
solved only using the same number of auxiliary restrictions, for 
example ( fTTT ). 

Below, we discuss two approaches to the extraction of & , & , 
pi, and di from A im . 



equal to the variance of a term Yjj' a U x jm 
duced applying restrictions 



44* = min 



in Eq. (|3), can be re- 



(15) 



on solutions ay of the system (0. 

As it was mentioned above, geometric field distortions (their 
variations in time), up to order k, are excluded along with image 
motion. 

This simple method was applied for the test reduction of 
four-hour series of frames obtained at VLT/FORS2 during one 
night (Paper I). By filtration of atmospheric image motion car- 
ried out with k up to 12, 300 pas precision of a single position 
for bright stars was achieved, which is near to the accuracy of 
photocentre measurements. Approximation Aj m = was valid 
due to a small residual DCR effect for frames obtained with the 
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Fig. 1. x, y residuals of Eqs. ([Pil l for 25 images computed 
with the assumption of "fixed" background star positions (left 
panel) in comparison to the exact solution (right) for epochs of 
reference and processed frames being separated by two years. 
Dot sizes are proportional to the stellar magnitude (only bright 
B <21.5 mag stars are shown). Note the difference in scales 
(1 mpx=100 pas). 



3.3. A case of fixed local reference star positions 

In the case of small background star displacements Aj m , and 
due to the averaging effect at summation, a term Yj/ a ij^-jm in 
Eq. ( fl~3T > can be discarded as negligibly small. In fact, this as- 
sumes a fixed position of local reference stars A jm = that leads 
to the most simple result A? - Xj m . For each target i, we are 
facing a system of 2M equations 

e + h x - h* = X im - e x - AM X 

e + W. - = Y im - e > - AM- (14) 
m= l,2...M, 

for determining ^ = |f - ef - AM* Q , f, = % - < Q - AM y m , p«, 
and di. This system can be solved by the least-squares method 
for each target object independently, with reference to its own 
subset Q;„, of reference stars, considering e m and AM; m as ran- 
dom (with respect to index m) errors with variances e? and cr 2 aV 
respectively. Errors e j m of reference star photocentre determina- 
tions add a further uncertainty to the phase of compution of X- m 
by means of Eq. A variance of this error cr^ = 'Z 1 j'af.s 2 ' m , 



3.4. Reduction for floating reference stars 

For the processing of sky images in the B-filter, the use of the 
technique descrived about is however limited, even for frames 
taken within the same night. In this case, an excess of posi- 
tional errors over s,- by a factor 3-5 results, because the under- 
lying assumption of A ]m - does not hold due to large (~2- 
10 mas/hour) chromatic displacement of images. For frames 
taken two years after the reference frame, which results in sig- 
nificant additional proper motion displacements, the solution of 
Eqs. (TBI produces unacceptably large residuals (FigQJ). 

The exact solution of Eqs. ( fTTITQl ) ensures a linearity and 
uniqueness of the transformation from X m to A m even for mov- 
ing (Ai m + 0) reference stars. Therefore, although X lm are deter- 
mined relative to local reference fields Q. im centered at the z'-th 
target and depend on the size Ri of these fields, the final solu- 
tions Ai m are consistent within the global reference field of size 
R and do not depend on the choice of initial Rj values. Since 
usually R > Rj, the transformation from X, m to A lm increases the 
atmospheric image motion component in /l„„, which should be 
compensated for by moving to larger k. 
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Note that due to conditions ([8j, both A, m and A- m + (fi, n are the 
solutions of a system ( [TBI if <£>,,„ are functions 



2 E wm f w 

w=l 



(16) 



with arbitrary coefficients E* m and E\ vm . Therefore, for splitting 
Ai m into components we can use 



£ x +h x -h x +tp x =A X +N X 

im /() ~ im tm tm 

m=\,2...M, i = l,2...N 



(17) 



where N; m = -e,„, - AA1 !m is a noise with variance sf + o~\ v 
Eqs. ( fTTT i can then be solved as a system of 2 2 N m conditional 
equations with AN + 2MN' unknowns gf, pi, dj, E* m , and 



^ wm • 



Because for most not very faint targets, cr at <sc e im , the least- 



squares solution is found with weights g,, 



„-2 



Also, due to 



the dominant contribution of bright stars, we can safely assume 
sin ~ lil&in wn i cn follows from the approximation given by 
Eq. (0 for these stars. 

Linear relations between some parameters of a system ( fTTI i 
do not permit to separate them completely. For example, system- 
atic variations of f?, p,, and d; quantities over a field imitate the 
behaviour of tp x m functions. Associations of this type are elimi- 
nated by applying N' orthogonalizing restrictions 



£ Pifmh = 0, Z diUh = 0, w = 1,2... AT' 



(18) 



The restrictions ( fT~8b permit us to recover only those components 
of pi, dj, and i* 1 ., that do not correlate with functions /„., . 

An essential drawback of a straightforward solution of Eqs. ( 
[T6l - fT8l ) is a large number of unknowns E x , m and E\ vm , that cause 
a certain instability of a solution. For example, the number of 
unknowns for Dec 2000 frames is 2MN' = 2 • 40 • 36 ~ 2900 
at k = 16, which is comparable to the number of measurements 
2 2JV m *9100. 

However, there is a simpler solution. Consider equation with 
a diagonal element E x vm of the normal system of Eqs. ( fTTI i: 



i€ox m 



(19) 



where summation is performed only over a limited sample a) m 
of stars measured in the m-th frame and, since FWHM»;const 
within the frame, the approximation g, m = 7, was used. Adding 
the sum 



A 



l m = J] {g+Kn-K&Mi 



(20) 



!£ nOt OJ,„ 



formed over stars not measured in the frame m to both sides 
of Eq. ( fT9l ), and taking into account ( fT~8b , we note that the first 
component in Eq. ( fl~9b turns to zero. Therefore 



(21) 



If a solution of the system ( fT3l for the m-th frame is found that 
fulfills the conditions 



Yj^mfvli^-Km W=\,2...N', 



(22) 



instead of those of Eq. ( fTTT ). the right part of the system given by 
Eq. ( |2T1 ) turns to zero, thus yielding a unique solution E wm = 0. 
Similar considerations for the y-component of the data yield 
E wm = 0, hence tp x n - = 0. The elimination of the large num- 
ber (2MN') of unknowns E x vm and E x m in Eq. ( fTTl l occurs due 
to imposing an equivalent number of restrictions d22l >. thereby 
increasing the effective degree of freedom of the system given 
by Eq. (TlTb to 2 2 N m - 4-N. Furthermore, this splits the system 
given by Eq. ( TPTb into A' independent subsystems of 2«, equa- 
tions (n, is number of observations of the i-th star) 



(23) 



which is easily solved for p,, and d, for each of the star ob- 
served independently. With this processing technique, A wm en- 
tering Eq. ( l22b are calculated by iteratively refining gf, , p, , and 
df. Iterations are however not required if each star is measured in 
all of the frames since in this case the identity A x vm = A\ vm = 
holds due to Eq. ( fTBl . 

Some comments should be made on the nature of £ values 
computed not in the absolute reference frame, but only relative to 
a limited set of nearby stars. Systematic differences between 
£ and values £ a t, s referring to the absolute or global frame, are 
given by 



fabs,i = 6 + 7 . G K f wi 



(24) 



which incorporates a polynomial function Yj W G w f w ; of N' or- 
der with coefficients G w . The values G w should confirm the re- 
strictions given by Eq. (TT~8b which now take a form 2;(^abs,i - 
YjwG w fm) 2 Ii = min, or Yti^h - mm - Eq. d24b thus is an ex- 
pansion of £ib s ,, into base functions and are the residuals 
of the expansion. The increase of N' (or k) leads to the increase 
of polynomial modes excluded from £ a t, s and thus to the loss 
of information retained in Reduction therefore should use the 
lowest k that still ensures acceptably small <r at . Above consider- 
ations, of course, are fully applicable to p, and dj parameters as 
to the remainders of corresponding polynomial expansions. 

4. Reduction of measured images 

An image obtained on 2 1 Dec 2000 at an hour angle near to the 
mean hour angle of other images was taken as a reference frame. 
Only stars in the central region of the sky with an R — 1000 px 
radius containing = 169 star images down to B = 24 mag were 
selected for the reduction. Peripheral stars which due to a strong 
jittering were often beyond frame boundaries thus were rejected. 
Nevertheless, the number of stars measured varied among the 
frames and reached down to N m » 90. 

We have taken advantage of the facts that the colours of the 
observed stars do not vary significantly on short time-scales (few 
days) and positions of most stars are only slightly affected by 
proper motions and parallaxes. Therefore frames were combined 
into series within annual epochs, which statistically improved 
the accuracy of model parameters determination. We therefore 
organized our data into two series of 39 and 25 frames for 
2000 and 2002-2003 epochs, respectively, each one represent- 
ing four nights. Astrometric reduction was carried out with dif- 
ferent choices of the parameter k, where the accuracy of results 
became sufficiently good at k — 16, which was adopted for the fi- 
nal reduction. Our iterative procedure involved the computation 
of gf, pi, as well as dj and the refinement of A,„, quantities 
calculated on the basis of the former parameters. The number 
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of outliers at each coordinate was limited to 1.25% of the total 
number of measured star images. 

The residuals u x m and u) m of the conditional equations 
Eq. d23l on x- and y- coordinates are uncorrelated (FigQJi) and, 
besides, substantially smaller than those obtained with the ap- 
proximation of "fixed" background star positions (FigQJ). We 
also did not find any correlation between m„„ and the stellar po- 
sition within the frame or the colour-dependent parameters p, 
and di. As a simple illustration for a further discussion, Fig[2] 
presents examples of u* residuals for stars of different magni- 
tude, from bright B = 18.5 to faint B = 23 mag. A difference in 
the data point scattering reflects the dependence of photocenter 
error on the light signal. No apparent correlation in time is seen. 
More extensive statistical analysis of M lm values is performed in 
Sect.5. 
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C m . A good consistency of C m and C' m values suggests that actual 
angles f m are really close to g' m . According to Jehin E. (private 
communication, 2006), this may occur when the LADC is not 
reset after finishing preceding observations of the same target. 
The final reduction was carried out with (1 + C,„)tan£ m rather 
than tankm- 
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Fig. 3. The LADC compensating displacement d; and atmo- 
spheric differential chromatic refraction p, for stars of B = 
18-24 mag (dots) with a linear (solid line) fit. The dot size 
is proportional to the stellar magnitude. Individual error bars are 
shown. 



Fig. 2. Example of u? residuals for 2000 epoch frames. 



4.1. Performance of the LADC 



Table 2. Relative corrections to tanf m : producing best accuracy 
of observations (C m ) and corresponding to a 10 mn delay in the 
LADC setting (C' m ) 



night 


Cm 


c m 


night 


Cm 


Cm 


2000 






2002-2003 






20 Dec 


-0.17 


-0.15 


29 Dec 


-0.06 


-0.11 


21 Dec 





-0.02 


30 Dec 


0.05 


0.09 


22 Dec 


0.26 


0.20 


31 Dec 


-0.11 


-0.15 


23 Dec 


-0.11 


-0.14 


02 Jan 


0.02 


-0.13 



Zenith distances £ m of the LADC setting used in Eqs. ( fT2l are 
not specified in the fits file headers. According to the technical 
documentation of the VLT, ( m is equal to the telescope zenith 
distance at the midpoint exposure time of the first science im- 
age of the series and is constant for all frames within the series. 
We however tried to use corrected values (1 + C m ) tan f m instead 
of tan (,„ in order to improve internal precision. Coefficients C,„ 
(Table |2]l that yield the best-achievable precision have been de- 
termined with an accuracy of about ±0.03 relative to some zero- 
point night, chosen as 2 1 Dec 2000. Closer examination revealed 
that (1 + C m )tan£,„ is approximately equal to tan^' n computed 
with the telescope zenith distance £' m 10 mn prior to the first ex- 
posure. Writing t,' m in the form tan^,', = (1 + C' m )tan£ m allowed 
us to introduce new values C' m and compute them in analogy to 



An excellent performance of the LADC is demonstrated by 
Figf3] which reproduces the distribution of p,- and d t (averaged 
over both epochs) with error bars based on the effective preci- 
sion of a single measurement (|26*] i. Blue stars are shown in the 
upper left corner. The distribution of stars is strongly concen- 
trated along a line p,- = - 1 2dj with a scatter that only slightly 
exceeds random errors. 

4.2. Effective precision of a single measurement 

Astrometric positions are affected by various sources of error. 
The most important are fundamental (photon noise in the image), 
of instrumental origin (optical field distortions, errors related to 
CCD, etc.), induced by the atmosphere (DCR effect, image mo- 
tion), or related to inefficient reduction techniques that should 
take into account and exclude largest errors. Without analyzing 
each source separately, some of which are rather efficiently fil- 
tered out at the phase of the reduction, we consider their total 
impact on the key system of Eqs. d23l ). whose solution yields fi- 
nal astrometric parameters £ A and £ v . This total error in Eqs. d23l 
was denoted as N m and includes e, m , AW(, m , and a sum of other 
residual errors not filtered by the reduction model. Ni„, includes 
both a random noise and components correlated in time. The 
presence of correlated components leads to strongly biased re- 
sults if the correlation time exceeds the length of a single series, 
while being shorter than duration of the observational compaign. 
In this case, processing a short series of frames will recover noise 
signature rather than the desired real astrometric signal. In next 
section we described results of a statistical analysis of N m based 
on measured residuals «,■„,. 

The variance crf m of the noise 7V„„ is the effective variance of 
a single measurement. Obviously, cr? m > sf m + o\ v The size of 
cr, m depends on the stellar magnitude the seeing #, m , the sky 
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Fig. 4. Effective precision of a single measurement cr (mean 
for x and y axes) for 2000 (circles) and 2002-2003 (triangles) 
epochs as a function of B. Centroiding error s at average seeing 
and background noise is shown by solid (first epoch) and dashed 
(second epoch) lines. A scale is 1 mpx=100 //as. 



level b m , the reduction parameters k and R, etc. In the following, 
the analysis simplifies by using the improved expression for e lm 
that, as compared to Eq. (0, more adequately takes into account 
variations of Q lm and b m , amounting to about +30%, and is valid 
in a wide (about 5 mag) range of stellar brightness. Considering 
asymptotic dependences of s on 9 and b given by Irwin ( 1985) 
for bright and faint images, we modified Eq. (|2]i to 



2.46 y/TiX 0o 



(25) 



The second component in parentheses is significant for faint 
B > 21.5 stars only. The Q m , or FWHM value, is available as a 
PSF model parameter (Sect.2). We have found that the variance 
of the normalized deviations «, m /e, m based on our new relation 
is almost independent from observation conditions. For that rea- 
son, the actual observed value of erf was obtained as a system- 
atic correction to e? : 



bn'^bn 



(26) 



where an average was taken over stars of approximately the 
same magnitude. Further averaging of erf with respect to i and 
m (over observations) yielded the effective variance of a single 
measurement cr 2 at average observational conditions during the 
considered epoch. Fig|4]presents both cr and the mean centroid- 
ing error s corresponding to average observing conditions as a 
function of B. Because cr matches well the errors s predicted by 
our numerical simulation, we conclude that the sum of all errors 
arising from different sources, including cr at , is small in compar- 
ison to e, so that, approximately, cr w s. This means that the 
astrometric precision depends only on errors of the photocentre 
determination which are random Gaussian and uncorrelated with 
time due their nature of arising from photon statistics. The latter 
property is very important since it forms the basis for substantial 
statistical improvements of the astrometric accuracy A(n,) for a 
sequence of n, frames with the increase of the number of mea- 
surements as 



where s is the photocenter error at average observation condi- 
tion. The estimate d27| i is approximate since assumes the best 
situation with a diagonal covariance matrix of the system ( f23l . 
The validity of Eq. d27| i is limited with respect to the maximum 
amount of frames n max , at which the value of A(n,) reaches a 
floor set by systematic errors. 

5. Testing the VLT temporal astrometric stability 

5. 1 . The Allan variance 

The simplest empiric validation of Eq. d27| > and illustration of 
how this law works at different is associated with use of 
Allan variance cr 2 KX of the residuals «,,„. This quantity (Pravdo 
& Shaklan ["l996l l is a function of time lag t < n,/2 expressed in 
frames 



1 



2(m + 1 - 2r) 



/j, + 1-2t 

z 



9=0 



p+q ■ 



r) 



(28) 



A(n,) = s/ yfni; n, < n n 



(27) 



ct^j(t) is a precision of the average of u m taken over t frames, 
and so corresponds to A 2 (t). For uncorrelated sets of M, m , the 
expected dependence cr 2 M ~ r _1 is consistent with Eq. ( l27l >. In 
effort to increase the range of time lags, we performed compu- 
tations of cr^j treating all epoch frames as a single series of 
length. For miscellaneous reasons, n,- varied from about QAM 
to 1.0M. Results for stars of different B magnitude classes are 
shown in Fig|5](two bottom graphs). Though, in general, plots 
for separate stars follow the expected law t _1 (straight lines with 
an ordinate s 2 at t = 1), it is difficult to make any conclusion at 
t > 10. To smooth statistical fluctuations and exclude the depen- 
dence on the star brightness, we normalized individual cr 2 Al by e 2 
(upper plot of Figj5]) and then performed averaging of results 
over all 169 stars. Obtained smooth functions (cr^/e 2 ) (black 
dots) indicate no evidence of large systematic error, and extend 
a validity of the approximation cr 2 M ~ r~' at least to t ~ 15. 
Hence Eq.d27l> is valid to about n max a; 15. The asymptotic floor 
of the Allan variance at large r is under the detection limit of this 
approach. We can claim however that with 15 frames, the result- 
ing precision of the average is about 50//as for bright stars. In 
Sect. 5. 3, we renew this discussion to derive more exact estimate 

of nmax. 



5.2. Distribution of residuals 

The null hypothesis Hq formally arising from above analysis is 
that the Ni„, component in Eqs. d23l represents an uncorrelated 
Gaussian noise. The ensuing consequences of this assumption 
are 

- the averaging law, Eq. ( l27l i. and 

- the absence of other noise components, including instrumen- 
tal errors, with a noticable magnitude. 

We applied several statistical tests of our hypothesis. In the 
first test, we formed a histogrammic distribution W of normal- 
ized deviations = V«/7( n i _ 2)ui m /cr im and compared it 
with the distribution Wo of these residuals expected when the 
data fulfill the hypothesis Hq. In order to obtain Wo, we per- 
formed some numerical simulations of observations with a typ- 
ical spacing of frames on hour angles, and adding a Gaussian 
noise with a unit variance in x and y to the right side of Eqs. (|23l . 
Observed distributions W were computed using all combinations 
of data sets available (2000 and 2003/2003 epoch frames, x and 
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Fig. 5. Allan variance as a function of time lag r for individ- 
ual stars (dashed lines) in comparison to the r _1 law (straight 
lines) in narrow magnitude ranges (two lower graphs) and the 
normalized variance ct 2 m Ie 2 for all stars (upper plot). The aver- 
age (cr^/e 2 ) taken over all stars is shown by black dots. The 
plots refer to x coordinate and 2000 epoch frames. 



y data, with faint, bright, or a whole sample of stars). FigJfJ, 
b present an example of such a histogram W for bright stars. 
Corresponding simulated distributions Wo are shown by smooth 
curves, for which the same binning as for the computation of W 
has been applied. Unlike W, m , the distributions Wo and W are 
not Gaussians and have a more strongly peaked shape. The con- 
sistency of observed and model distributions by visual inspec- 
tion is rather good. Numerically, it was estimated by \ 2 values 
computed from differences between W and Wo in fixed binned 
intervals containing at least 5-10 data points. The^ 2 values with 
corresponding degrees of freedom (Dof) are listed in Table[3]for 
faint, bright, or all stars, respectively. The^ 2 -test based on these 
data does not reveal a deviation of W from Wo at significance 
levels below 20%, which does not allow to reject our hypothesis 
Hq. Due to a limited number of data points, this test is however 
only sensitive to a central ±2cr region of the distribution. 



Table 3. x 2 values and Dof characterizing differences between 
W and Wo 





faint bright all stars 




2000 epoch 


x 2 > x ly 

Dof 


20.1 /21.5 23.0/22.7 24.4/28.6 
22 22 24 




2002/2003 epoch 


Dof 


20.8/11.6 25.5/17.2 27.8/27.1 
18 20 22 
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Fig. 6. Histograms of individual normalized residuals u', m in x 
(a) and in y ( b) for actual (steps) and simulated (smooth lines) 
observations, on 2000 (solid) and 2002 epochs (dashed lines); 
the same for normalized night average residuals Ui P /o- p (c and 
d). Plots are drawn for bright B < 20.5 stars. 
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Fig. 7. The diffences D„ between the observed and model cu- 
mulative frequency distributions of normalized residuals u' lm for 
2000 epoch measurements in y (bright stars) 

We applied the more powerful Kolmogorov-Smirnov test to 
detect a potential difference between the observed and simu- 
lated distributions. For this purpose, we computed cumulative 
frequency distributions of observed and model normalized resid- 
uals u'im, and then obtained the difference D„ between two dis- 
tributions. A case with largest deviations for 2000 epoch mea- 
surements in y with 1721 data points (bright stars only) is shown 
in Fig|7]where the differences D n are plotted as a function of the 
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normalized deviation. The Kolmogorov-Smirnov test based on 
this plot (in this worst case max|D„| = 0.024) shows no signif- 
icant difference between the model and measured distributions 
at 20% confidence level. Both tests applied thus do not reject 
hypothesis Hq of a really normal distribution of W,,,, at approxi- 
mately the same significance levels. 

For another test, we considered nightly averages Ui p = 
Y iij m crT 2 1 V o"7 2 involving all observations of the 2-th star 

mep mep 

within a single p-th night. Fig|6]c, d show observed and model 
frequency distribution of normalized residuals Ui p /cr p , where cr p 
is the error of the determination of Ui P . A good agreement of the 
histograms of the observed statistics (both for bright and faint 
stars) with their expectations from model simulations, a symme- 
try, and the absence of large deviations (normally less than 1% of 
Uip/cTp values deviate at about 3 units) demonstrate excellent sta- 
bility of the VLT instrumental system during time spans of a few 
nights. Compatability of the observed and model distributions is 
also confimed by the Kolmogorov-Smirnov- and^ 2 -tests. 

5.3. Integral statistics 



tides are related by S^(obs) = (^ 2 of /Dof)S 2 (mod) + S 2 ys . For 
statistical reasons, it is probable that the sampled value 5^(obs) 
is equal or below of its mathematical expectation S 2 (mod) + S 2 ys . 

The largest value of S 2 ys , according to Table |4] can exist in the 
x-axis measurements of bright stars (2000 epoch), for which we 
found S 2 ys /S A (mod) < 0.19 at 20% confidence level. This is 
a ratio of a systematic to random error component for a single 
night series represented, in average, by 8 frames. In other words, 
the bias related to systematic errors is equal to a random com- 
ponent expected in positions averaged over 8/0.19 = 40 frames. 
This is an estimate of n max introduced in Eq.((27j, more precise, 
as compared to that obtained in Sect.5.1. For y-coordinate, we 
obtain n max = 65. Frames of 2002 yield n max = 120 and 50 for 
the x and y axes correspondingly. With n max = 50 adopted as a 
reliable limit, precision of a night series is about 30 ^as. 

5.4. Correlations 

A final important test was carried out to investigate correlations 
in JV„„ which can negatively affect the averaging law d27l ). For 
that purpose, we computed the autocorrelation function of nor- 
malized residuals 



Table 4. Integral statistics So and Sa for bright B < 20.5, faint 
B > 20.5 stars, and for simulated observations 



r(f) = H(f){(u im Jcri mi )(ui m ,/<Ti m ,)) 



(29) 





S£ S y n error 


S A S y A error 




2000 epoch 


model 
bright 
faint 


0.997 0.985 

0.992 0.989 ±0.021 

0.982 0.988 ±0.019 


0.749 0.839 

0.784 0.849 ±0.064 

0.761 0.787 ±0.054 




2002-2003 epoch 


model 
bright 
faint 


0.994 0.973 

0.973 0.961 ±0.029 

0.989 0.952 ±0.026 


0.738 0.845 

0.722 0.872 ±0.067 

0.744 0.874 ±0.055 



We also considered integral statistics St 



1 1 y («im-»ip) \ 

\„,-4^ o- 2 



where averages are taken over all stars and nights, and S 2 A = 

(ff£S? /cr 2 ), where an average is taken over all stars. These 

p 

two statistics are either sensitive to the variation of systematic 
errors within a given night (So) or from one night to another 
(Sa), and therefore are a powerful indicator of their presence. 
The quantities So and Sa, calculated separately for both coordi- 
nates x and y, are listed in Table [4] for each epoch of observa- 
tions together with the respective expected values 5'o(mod) and 
SA(mod), obtained from numerical simulations where random 
Gaussian noise W, m has been assumed. Errors due to a limited 
number of measurements indicate an 80% confidence interval 
for possible deviations of observed values from our model as- 
sumption. Observed statistics match well the model within the 
error limits. 

The fact that we do not find an excess in the measured So 
and Sa values means that no or very small extra noise except 
Ni m could be present in measurements. It is useful to estimate 
its upper limit. Suppose that the observed positions, besides of 
Ni m , are affected by some small systematic error that is constant 
within a night but varies between nights with an amplitude pecu- 
liar to a certain star. The observed value of S^obs) then exceeds 
its model expectation S 2 A (mod) by a small amount S 2 ys which is 
the average signal variance for a given star sample. These quan- 



of the argument r = ;«2 - m\, which is equal to the difference 
of frame indexes and used as a time lag. Here, H(t) — 1 - t/M 
is the weight function that compensates the increase of statisti- 
cal variations at large r. Averaging is performed over the images 
of all stars. A corresponding computation of the expected cor- 
relation function r(r) was carried out with the model observa- 
tions containing white noise with a unit variance. The observed 
r(j) values computed for both coordinates x and y for bright 
B < 20.5 mag stars are shown in Fig[8] along with the model 
function r(r) (mean for x and y) and 1-cr limits for statistical 
scatter of individual deviations. Oscillations in r(f) plot inherit 
a compound four-night structure of the series. Observed data in 
general follow the expected dependence, though with few iso- 
lated deviations over 2 sigma. 
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Fig. 8. Observed autocorrelation function in x (circles) and in y 
(squares) as compared to the model function (solid line) with a 
1-cr tolerance limits (two dotted lines) and approximation OOb 
(dashed line) 



Both r and r are not identically zero as it could be expected 
for uncorrelated measurements. The observed systematic neg- 
ative bias of correlations, however, is not due to instrumental 
errors but has a theoretical background (Lazorenko l 19971 ). In or- 
der to see this, consider Eqs. d23l for the i-th star. When the 



10 



P.F.Lazorenko et al.: High-precision astrometry on the VLT/FORS1 at time scales of few days 



number of frames M is large and the distribution of frames with 
hour angle is sufficiently random, all parameters in the system 
given by Eq. ( 1231 become uncorrected. Then is approxi- 
mately equal to the average of (A? - h x im + h* Q + Mm) taken 
with respect to the index m. In this approximation, one finds 

u im ~ (^L ~ Kn + h m + N '""> ~ €f- The term 8 iven in parentheses 
formally represents a solid time series of the "time" argument 
m with m = 1 . . .M measurements containing noise. With this 
approach, the u x , values are the remainders of a time series af- 
ter subtraction of the fitting polynomial (namely a constant g x in 
our case). The bias of the autocorrelation function shape due to 
subtraction of the best-fitting polynomial for measurements lim- 
ited in time was also previously studied by Lazorenko ( 11997b . 
The autocorrelation function r' of series remainders takes a par- 
ticularly simple form if the power spectrum of measurement er- 
rors resembles a white noise (which is the case) with a unit vari- 
ance and constant rate of data sampling. In this case, one obtains 
(Bakhonskiet al. [T997i l 

2c 

r'(r) = --^-sinc(27rcT/M) (30) 

where g = 0.56, if the subtracted polynomial is presented by a 
constant, and sinc(z) = sin(z)/z. One can see that r' « for large 
M only (long series) while for small M (limited number of sam- 
pled data) its value is absolutely large and negative. Due to sim- 
plification of our considerations, the measured and model func- 
tions are not exactly fitted by r', but good enough however to ex- 
plain the statistical origin of the observed negative bias (Fig. [8]). 
We conclude that no correlation is present in TV,,,, at scales of 
several frames. 

Table 5. Observed r and model r mo d normalized correlations. 



t' 


r v 


'"mod 


-V -V 

mm 


1-sigma error 




2000 epoch 


1 


-0.17 


-0.18 


-0.22 -0.24 


±0.05 


2 


-0.17 


-0.10 


-0.20 -0.24 


±0.07 


3 


-0.42 


-0.38 


-0.25 -0.23 


±0.09 




2002-2003 epoch 


1 


-0.15 


-0.14 


-0.26 -0.25 


±0.06 


2 


-0.19 


-0.20 


-0.32 -0.23 


±0.07 


3 


-0.35 


-0.26 


-0.13 -0.19 


±0.10 



Correlation at time scales of few days were studied using 
normalized nightly-average residuals u ip /cr p . We computed cor- 
relations 

f(j') = ((u ip Jo- pi )(u ipi /o- p2 )) (31) 

as a function of the time lag t' = pi - pi between two nights 
with indexes p\ and pi, expressed in days. The values of f for 
the observed data and corresponding model value r mo d are given 
in Table[5] Apparently, there is a strong negative bias in f due to 
a very limited number of data points in our 4-night series. This 
bias, however, is well-modelled using uncorrelated noise N lm , 
which confirms previous conclusions that there are no traces of 
instrumental signature in VLT observations at time scales of few 
nights, and supports the validity of the averaging law (l27l >. 

5.5. Two-year stability 

Given only two epoch measurements, no conclusions on a long- 
term astrometric stability of the VLT can be made based on 
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Fig. 9. Difference in the chromatic parameters p, and dj between 
the two epochs, normalized by the error of this difference com- 
putation (dots). Ellipses are drawn at 1% (dashed line) and 5% 
(solid) significance level. A straight line shows the direction of 
point displacements due to the actual change of star colours. 

positional information. Considering extreme importance of this 
problem, we present here indirect analysis of this issue based, in- 
stead of positions, on a comparison of chromatic p, and d, param- 
eters computed from independent runs of 2000 and 2002/2003 
epochs. If one assumes constant colours for the observed stars, 
these parameters should coincide within error bars. A change 
of the VLT characteristics, depending on stellar colours, should 
result in an additional difference. Fig|9]shows normalized differ- 
ences between p, and dj between the two epochs, where ellipses 
refer to significance levels of 1% or 5%, respectively. Some 
points represent outliers, testifying the presence of systematic 
errors. The largest residuals for bright stars in the right bottom 
however occur along a line p,- = - 1 .2af,- that is peculiar to the dis- 
tribution of chromatic parameters shown in Fig. [3] It is therefore 
very likely that these deviations are caused by a change in stellar 
colour. Some other large deviations appear to be related to stars 
in outer regions of the observed field where the method seems to 
give biased results. Except for these outliers, the distribution of 
points in Fig. |9]resembles a Gaussian, which is an indicator of 
a rather good astrometric stability of the VLT even at very long 
time scales. It should be stressed that above test is insensitive to 
non-chromatic type of errors and therefore is only indicative. 

6. Conclusion 

Our results reveal an exceptionally good astrometric perfor- 
mance of the VLT and its camera FORS 1 . Thus, the precision s 
of the position of a bright star reaches 200 ^as for a single mea- 
surement, which fairly well corresponds to the former 300 pas 
estimate obtained for FORS2 camera (Paper I). Both cameras 
thus give similarly good results. The term "bright star" refers to 
unsaturated images containing 1-3T0 6 electrons and may cor- 
respond to different stellar magnitudes depending on the expo- 
sure, filter, and seeing. Most importantly, we found negligibly 
small systematic errors of instrumental and other origin. In fact, 
no traces of these 'dangerous' errors were found at time intervals 
up to 4 days. Due to this fact, the precision for a frame series im- 
proves as e/ -\fh~i, and at m = 50 reaches 30 pas for FORS 1 and 
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40 yuas for FORS2. With use of FORS2, the observation time 
needed to obtain this accuracy is about 1 hour if short 15 s expo- 
sures and a 2 x 2 binning of pixel reading is used. We conclude 
that the VLT with cameras FORS1/2, due to its enormous col- 
lecting light power, fine optical performance, and effective aver- 
aging of wave-front distortions over a large aperture, is a pow- 
erful instrument that can be used efficiently for high-precision 
astrometric observations of short-term events, in particular, of 
planetary microlensing. 
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